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Abstract 



In this formahsm, various logarithmic corrections are organized to give the 
Wilson evolution from the W boson mass down to the characteristic scale 
of a decay process, which is of order of the b quark mass, and the Sudakov 
evolution from the characteristic scale to a lower factorization scale of order 
^QCD- With large logarithms organized, the b quark decay amplitudes are 
evaluated reliably in perturbation theory. Below the factorization scale, QCD 
dynamics is regarded as being nonperturbative, and absorbed into meson 
wave functions. Because of their universality, the heavy-meson wave func- 
tions determined in this work, can be employed to make predictions of other 
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decay modes. We also observe that the dependence of heavy meson wave 
functions on intrinsic parton transverse momenta plays an important role in 
the explanation of data. 



I. INTRODUCTION 

Recently, a PQCD factorization theorem for exclusive heavy-meson decays has been 
developed by one of us [0], and applied to the nonleptonic B -^ D^*^n(p) decays and 
the penguin-induced radiative B — > K*'y decays 0. These decays involve three scales: the 
W boson mass M^, at which the matching conditions of the effective weak Hamiltonian to 
the full Hamiltonian are defined, the typical scale t of a hard amplitude, which reflects the 
specific dynamics of a heavy-meson decay, and the factorization scale 1/6, with b being the 
conjugate variable of parton transverse momenta. The dynamics below 1/6 is regarded as 
being completely nonperturbative, and parametrized into meson wave functions (f){x,b), x 
being the momentum fraction. Above the factorization scale, PQCD is reliable and radiative 
corrections produce two types of large logarithms: hi{My//t) and ln(t6). The former are 
summed by renormalization-group (RG) equations to give the evolution from My/ down to 
t described by the Wilson coefficients c(t), while the latter are summed to give the evolution 
from t to 1/6. 

There also exist double logarithms In (P6) from the overlap of collinear and soft diver- 
gences in radiative corrections to meson wave functions, P being the dominant light-cone 
component of a meson momentum. The resummation of these double logarithms leads to 
a Sudakov form factor exp[— s(P, 6)], which suppresses the long-distance contributions in 
the large 6 region, and vanishes as 6 = 1/A, A = Aqcd being the QCD scale. This factor 
guarantees the applicability of PQCD around the energy scale of the 6 quark mass, which 
is our motivation to develop a PQCD formalism for heavy-meson decays. For a detailed 
derivation of the relevant Sudakov form factors, we refer the readers to |Q,|^ . With all the 
large logarithms organized, the remaining finite contributions are absorbed into the hard 6 
quark decay amplitude H{t). 

We do not claim that contributions to B meson decays are purely perturbative, but that 
perturbative contributions are important, and it makes sense to compare PQCD predictions 
with experimental data. It has been shown that about 70% of the full contribution to the 



decay B -^ Dn arises from the region with a^ < 1 (os/vr < 0.3) in the above PQCD formal- 
ism . Since energy releases involved in B meson decays are located in the transition region 
from nonperturbative to perturbative QCD according to the analysis in 0, nonperturbative 
approaches, such as QCD sum rules |0], also give sizeable contributions. In fact, B meson 
decays are a subject on which perturbative and nonperturbative approaches complement 
each other. 

A three-scale factorization formula for exclusive B meson decays possesses the typical 
expression, 



c(t) H{t) (g) 0(x, b) (g) exp 



-siP,b)-2f ^7(«.(/^)) 

Jl/b U 
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where the exponential involving the quark anomalous dimension 7 = — a^/vr describes the 
evolution from t to 1/6 mentioned above. All the convolution factors in the above formula, 
except for the wave functions 0(x, 6), are calculable in perturbation theory. The wave 
functions, though not calculable, are universal, since they absorb the long-distance dynamics, 
which is insensitive to the specific decay of the b quark into light quarks with large energy 
release. The universality of nonperturbative wave functions is the fundamental concept of 
PQCD factorization theorems. It can be proved rigorously, but the proof is not the point of 
this work, and will not be supplied here. 

Because of the universality, the strategy of PQCD factorization theorems is as follows: 
evaluate all perturbative factors for some decay modes, and adjust the wave functions such 
that predictions from the corresponding factorization formulas match experimental data. At 
this stage, the nonperturbative wave functions are determined up to the twists and orders of 
the coupling constant, at which the factorization formulas are constructed. Then evaluate 
all the perturbative factors for another decay mode. Input the extracted wave functions 
into the factorization formulas of the same twist and orders, and make predictions. With 
this strategy, PQCD factorization theorems are model independent and possess a predictive 
power. 

For light-meson wave functions, such as those of vr, p, and K^*^ mesons, various ex- 



perimental data and QCD sum rules p|-|TO| have provided reliable information. For heavy 
mesons, such as B and D^*\ the knowledge is still poor. In this paper we shall determine the 
heavy meson wave functions from the available experimental data. Since two-body nonlep- 
tonic decays involve the maximal energy release that guarantees the applicability of PQCD, 
and their data are abundant, the decays B -^ D^*^7i are ideal processes for extracting the 
B and D^*^ meson wave functions. However, both of the B and D^*^ wave functions are 
unknown, and thus we need the additional data of the B —* K*'y decay |ll| in order to 



determine the B meson wave function first. After determining the B meson wave function, 
we apply it to the analysis of the B -^ D^*^n decays and obtain the D and D* meson wave 
functions. 

The factorization formulas for the B -^ K*'j and B -^ D^*^n decays are presented in 
Sec. II and Sec. Ill, respectively. The numerical analysis is performed in Sec. IV. Section V 
is the conclusion. 

II. THE B -^ K*j DECAY 

The effective Hamiltonian for the flavor-changing b ^ s transition is given by [|1^] 

i/eff = -^V*Vt,j:Q{^^)0,{fi), (2.1) 






with 



Oi = {SiCj)v-A{Cjbi)v-A , O2 = {SiCi)v-AiCjbj)v-A , 

O3 = iSibi)v-A J2(^J^j)v-A , O4 = {Sibj)v-A J2(^J^i)v-A , 
q q 

O5 = {Sibi)v^Aj2(^J^j)v+A , Oq = {Sibj)v~A^iqjqi)v+A , 



e 
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O7 = ^^Sia'"'{msPL + rrihPRJbiF^, , 



Os = -£^rs,a^'im,PL + m,PR)T^^,G%, (2.2) 

i, j being the color indices. It has been observed that the contributions to the B -^ K*'y 
decay rate come mainly from the operators O2, Oj and 08 0.02 contributes via a charm- 



quark loop, from which the photon and an off-shell gluon are radiated. This off-shell gluon 
is then reabsorbed by the spectator quark in the B or K* meson. This b -^ sg*^ vertex has 
been computed in [P,p!^. 

The momenta of the B and K* mesons in light-cone coordinates are written as Pi = 
(Ms/V2)(1,1,0t) and P2 = (MB/y2)(l, r^, Ot), respectively, with r = Mk^/Mb- The B 
meson is at rest with the above parametrization of momenta. We define the momenta of light 
valence quarks in the B and K* mesons as ki and ^2, respectively, ki has a minus component 
ki , giving the momentum fraction Xi = ki / P{ ^ and small transverse components ki^. k2 
has a large plus component /cg', giving X2 = k2 / P2 ■, and small k2T- The photon momentum 
is then P3 = Pi — P2, whose nonvanishing component is only P^ . 

The B -^ K*'j decay amplitude can be decomposed as 

M = e; ■ e*^,M' + it^,+.e*^^e*^M'' , (2.3) 

with e^ and e^-* the polarization vectors of the photon and of the K* meson, respectively. 
The total rate of the B — > K*'^ decay is given by 

" ^"'" (|M^P + |M^n . (2.4) 



SttMb 
We further decompose M^ and M^ as 

M' = Mi + Mij + Mi, (2.5) 

where i = S or P, and the terms on the right-hand side represent the contributions from 
the operators O2, O7, and Og, respectively. 

The Sudakov resummation of the large logarithmic corrections to the B meson wave 
function 0b and to the K* meson wave function (px* leads to the exponentials exp(— ^b) 
and exp{—SK*), respectively, with the exponents 

Ssit) = s{xiP,-, hi) + 2 f %{asm)) , 
SK'it) = s{x2P^M) + s{{l-X2)P^M) + '^ f —liasm)) . (2.6) 

Jl/b2 U 



The variable 61 (62), conjugate to the parton transverse momentum kiT {k2T), represents 



the transverse extent of the B (K*) meson. The exponent s is |p!^ , p!5 
s{Q,b) 



i/b /i 



lnpM(a,(/i)) + 5(a,(/i)) 



(2.7) 



where the anomalous dimensions A to two loops and B to one loop are 
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A = Cp^ + 

TT 



3 TT V 2 



67 7r2 10 
Y~y ~ 27 



2 /e'^^\ 



(2.8) 



with Cp = 4/3 a color factor and 7^; the Euler constant. The one- loop expression of the 
running coupling constant, 



asifJ,) 



(2.9) 



/?oln(/iVA2) 

is substituted into Eq. ( |2.7| ) with the coefficient Po = (33 — 2nf)/3. The Sudakov factor in 
Eq. ( |2.6| ) suppresses long-distance contributions from the large b region, and improves the 
applicability of PQCD to B meson decays. 

We quote the factorization formulas for Mf [^] in terms of the overall factor 

Gpe. 



r{o) 



V2 



TC 



V*VtbCFM% . 



(2.10) 



The amplitudes contributed by O2 are written as 

A rl rl — x rl /"l/^ 

j^s'_-p(o)_ (ijr dy dxidx2 / 6(i60^(xi, 6)0x*(a^2) 

3 JO Jo Jo Jo 

Xas(t2)c2(t2) exp[-S'B(t2) - SK*it2)]bi=b2=b 

X [(1 — r^ + 2rx2 + 2xi)y — {rx2 + 3xi)(l — x)] 

,, (l-r)(l-r^)x2X f—^ 

^a;y(l-r>2M|-M,2^^^^^'V'^^'^^ ' 

M2 =r^°-'- dx dy dxidx2 / bdb(f)B{xi,b)(l)K*{x2) 

3 JO JO JO JO 

Xa^(t2)c2(t2)exp[-5'B(t2) - 5'i^*(t2)]6i=62=6 



f2.111 



((1 - r)(l - r^) + 2r^X2 + 2xi) y 

— (r(l + r)x2 + (3 — r)xi) (1 — x)] 

(1 — r'^)x2X 
^ xy{l - r^)x2Ml - Mf 



H2{Ab,J\Bl\b), (2.12) 



with 



A^ = X1X2MI , 



B 



B'^ = xiX2Mi, - -^^(1 - r')x2Mi, + 



\ — X x{l — x) 



t2 = max(A,^|5||,l/6), (2.13) 

and Mc being the charm quark mass. ^2 is the characteristic scale of the hard amphtude 



H2{Ah, J\Bl\b) = Ko{Ab) - Ko{J\Bl\b) , B^ > , 



TT 



KoiAb) - i'-Hl^\^\h) , Bl<Q. (2.14) 

The amphtudes from 0^ with a hard gluon attaching the quarks in the B meson or in 
the K* meson are 

M^ = -Mf = r(°)2 / dxidx2 / bidbib2db2(l)B{xiM)(t>K*{x2) 
Jo Jo 

xasih)C7itr)exp[-SBih) - SK^h)] 
x(l - r')LH!r''\Ab2, ^/Wlbi, ^/Wlh) 

+ [l + r+{l-2r)x2]H!,'\Ab,,CM,Cjb2)\ , (2.15) 

with 

Bj = {xi - r'^)m\ , Cj = X2m\ , 



tr = max(A, V|5?|, C7, 1/61, 1/62) . (2.16) 

The hard functions are given by 



r(a). 



H!,''>{Ab2,Jmbi,Jmb2 



K,{Ab2)h{J\B^,\b,, J\B^,\b2) , B',>0, 



K,iAb2)h'iJmb,, J\B^b2) , 5? < , (2.17) 



with 



h-- 


= e{b, 


- h)Ko{^mh)Io{^\B^j\h) + (61 ^ 62) , 


h'-- 




e{h - 62)iJi'Vl^7 h)M^Jmb2) + {h ^ 62) 



f2.18) 



and 



Hi''\Ah, Crh, Crb2) = K,{Ab,)h{C,b^, Cjb,) . 



(2.19) 



The relation Mj = —Mf reflects the equality of the parity-conserving and parity- violating 
contributions induced by O7. 

The amplitudes from Os, where the photon is radiated by each quark in the B and 
mesons, are 



K' 



[^2 



1 /"i /"iM 

Ma =—T^^'— / dxidx2 / bidbib2db24>B{xi^bi)(l)K-'^ 
3 JO Jo 

xas{h)Cs{h) ex^[-SB{h) - SK*{h)] 
xUl-r'' + xi){rx2 + xi)H^^\Ab2,B^b^,Bsb2) 
+ [(2 — 3r)a;2 — Xi + r(l — X2){rx2 — 2rxi + 3xi)] 

xHf\Ab,,C^b,,C^b2) 
+ (1 + r)(l - r2)[(l + r)xi - rx2]H't\^\A^\b2, Dsbi, Dsb2 
-[(1 - r2)((l - r^){2 - X2) + (1 + 3r)(2a:2 - Xi)) 



+2r^X2{x2 - x,)]Hi''\,/\A^\bi, Esbi, E^b2) \ , 

M[ = T^^'- dxidx2 bidbib2db2(j)B{xi,bi)(j)K*{x2) 

3 JO JO 

xas{ts)C^{ts) eM-Ssits) - SK^it^)] 
A{l-r'' + xi){rx2 + xi)Ht\Ab2,B^bi,B^b2) 
+ [(2 — 3r)a;2 — Xi — r(l — a;2)(ra;2 — 2rxi + 3a;i)] 



xHt\Ab^,CsbuC^b2) 
+ (1 - r)(l - r2)[(l + r)xi - rx2Wi\^J\A^\b2, Dsb,, D^b2] 
-[(1 - r2)((l - r2)(2 + 0:2) - (1 - 3r)xi) 



(2.20) 



rW, 



-2r'x2{x2-x{)]Hl^\J\A'^\h,E^huE^h2) , (2.21) 



with 



A'^ = (1 - r^){xi - X2)Ml , Bl = (1 - r^ + xi)M| , 

Cl = (1 - X2)M| , D^ = (1 - r^)xiMl , ^2 = (1 - r^)x2Ml , 



tg = max(A, ^\A% Bs, Cs, Ds, Es, 1/h, I/62) • (2.22) 

The hard functions are given by 

Hlt\Ab2,Bsh,Bsb2) = Hi'\Ab2,Bsh,Bsb2) , 
Hi'\Abi,Csbi,Csb2) = Ko{Ab,)h'{Csbi,Csb2) 



Hi'\J\A^\b2,Dsb,,Dsb2 



Ko{^\A'^\b2)h{Dsb^,Dsb2), A">0, 
i^H!,'\^\A^\b2)h{Dsb,,Dsb2) , A"<0 



rW, 



Hr{J\A'^\bi,Esbi,Esb2 



Ko{^\A'^\bi)h\Esbi,Esb2), A'^>0 



t-H','>{^\A'^\b,)h\Esb,, Esb2) , A"<0. (2.23) 

Compared with the factorization formulas presented in 0, we have included the intrinsic 



(nonperturbative) b dependence [0 of the B meson wave function (f)B ■ The b dependence 
appearing in the Sudakov factor is regarded as being perturbative. The intrinsic b depen- 
dence, providing additional suppression in the large b region, is more essential for the B 
meson, in which soft dynamics dominates and Sudakov suppression is weaker. This depen- 
dence for the fast recoiling K* meson is less important and neglected. In fact, the intrinsic 
b dependence of the K* meson wave function is not yet fully understood. For an analysis 
of a single process such as the decay B -^ K*'y p, the effect of the intrinsic b dependence 
is not significant, since it can always be compensated by adjusting the x dependence of (pB 
under the current accuracy of experimental data. However, this dependence turns out to be 
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essential, when one intends to explain data of several B meson decay modes simultaneously. 
This observation will become clear in Sec. IV. 

Another modification is that we have chosen a single hard scale t for all diagrams con- 
tributing to an amplitude, which is defined as the maximal scale appearing in these diagrams. 
In [^ t is chosen as the hard gluon momentum of each diagram, and thus differs among var- 
ious diagrams contributing to an amplitude. This modification simplifies the factorization 
formulas for the decay B -^ K*'y. The difference between these two schemes arises at the 
level of next-to-next-to-leading logarithms, if B meson decays are indeed dominated by hard 
gluon exchanges. 

III. THE B -^ D^*h DECAYS 

In this section we quote the factorization formulas for the nonleptonic B -^ D^*^tx decays, 
whose effective Hamiltonian is given by 

^eff = '^v:,v,,[CMO^{^^) + CM02{^^)] , (3.1) 

with the four-fermion operators 

Oi = idiUi)Y^A{cjbj)v-A , O2 = {diUj)v-A{cjbi)v-A ■ (3.2) 



Note the interchange of the definitions of Oi and O2 compared with those in Eq. ( p.ll) . 
Similarly, we assign the momenta of the B and D^*^ mesons as Pi = {Mb/V^){1, 1, Ot) and 
P2 = (MB/v^)(l,r^,OT), respectively, with r = M^m/Mb, M^h*) being the mass of the 
D^*^ meson, e* is the polarization vector of the D* meson. The momentum fractions Xi and 
X2, and the transverse momenta ki^ and \<.2t are defined similarly as in Sec. II. The pion 
momentum P3 = Pi — P2 is the same as the photon momentum in Sec. II. One of its valence 
quark carries the fractional momentum x^P^^ , and small transverse momenta ka^. 

In the present analysis we regard the D^*^ meson as a heavy meson. Hence, the Sudakov 
form factor exp(— S'^jc.)) from the resummation of the double logarithmic corrections to the 
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D^*^ meson wave function (j)jj{*) has the same expression as exp(— S'b) but with different 
kinematical variables. The Sudakov form factor exp(— 6*7^) associated with the pion wave 



function (pT^ contains the exponent 17 



S^it) = sixsP.i, 63) + s((l - xs)P^, 63) + 2 r ^7(as(/i)) • (3.3) 

The decay rates of -B -^ D^*^n have the expression 

^ :GiiK.nK.rM|ii^^iMr, (3.4) 



1287r ^' ' ' ' ^ r 
where i = 1, 2, 3, and 4 denote the modes fi" -^ D^-r-, B^ -^ D+vr", R- -^ D*^7t-, and 
i?'^ -^ D*+7r^, respectively. The decay amplitudes Aii are written as 

Ml = A[(l + r)^+ - (1 - r)^_] + /^U + -Mext + -Mint , (3.5) 

M2 = /.[(I + r)e+ - (1 - r)^-] + Mexc + A^cxt + A^exc , (3.6) 

Ms = ^M(l + r)a - (1 - r)irU + U)] 

+fD4L + M:^, + Ml, , (3.7) 

M, = ^/.[(l + r)U - (1 - 0(^U + ^s)] 

+/bCxc + -M:xt + -m:., , (3.8) 

where /s, /^(.), and /^ are the -B meson, D^*^ meson, and pion decay constants, respec- 
tively. The form factors C,i, i = +, — , V, Ai, A2, and ^3, denote the factorizable external 
PF-emission contributions. The form factors .^j^^ and ^^*), denote the factorizable inter- 
nal ly-emission and PF-exchange contributions, respectively. The amplitudes A^ixt) Ml^L 
and A^ixl; represent the nonfactorizable external ly-emission, internal ly-emission, and W- 
exchange contributions, respectively. 

The form factors ^i, i = +, Ai and ^3, and $,j, j = — and A2, are given by 

/•l /"OO 

^i = IGirCFVrMl dxidx2 bidbib2db2(f)B{xi,bi)(j)j;)(,){x2,b2) 
Jo Jo 

xas(t)ai(t) exp[-SB{t) - 5'£,{,)(t)] 

X [(1 + (iX2r)h{xi,X2, bi, 62, m) + {r + (lxi)h{x2, Xi, 62, ^i)] , (3.9) 

12 



IGtvCfVtM^ dxidx2 bidbib2db2(f)B{xi,bi)(j)j;)(,){x2,b2) 
Jo Jo 

xasit)ai(t) exp[-SB{t) - Sj;)(,)(t)] 

X [CjX2rh{xi, X2, bi, 62) + CjXih{x2, xi, 62, &i)] , 



(3.10) 



with the constants ITS 



c+ = c; 



' 2 V^+1 



77-1 / 1 

^"2 



C- = -c 



1 7/+1 / 3 

^-2 + V^Fl'v^"2 



a, 



2 — T] — \/rj'^ — 1 



2(7/ + 1) 



Cai 



2(77 + 1) 



Ca. = , C 



V 



A2 



c 



A3 



77-2 






(3.11) 



2 2v/^^^^T ' ^^^^ 2V^^^^T ' 

where the velocity transfer 7/ takes the maximal value 7/max = (1 + 7'^)/(2r) here. The form 
factors Q*l and ,^^*). are written as 

^i„t = IQ-kCeVtMI / (ixifixs / bidbibsdbs(f)B{xi,bi)(j)^{x3) 
Jo Jo 

Xas(tint)a2(tmt)exp[-S'ij(tint) - ^^(imt)] 
X (1 + X3(l - r2))/lint(Xi, X3, 61, 63, 77iint) 

+Cint a;ir^/iint(a;3, xi, 63, 61, mint) , 

("1 rco 

Qll. = IGnCpVrMl dX2dX3 62rf&2&3C?&30D(*)(a;2, &2)07r 

JO Jo 

xas(texc)a2(W)exp[-5'£,{,)(tcxc) - 5*^(^0x0)] 

X [(^3(1 - r^) - Ccxc'"^)^cxc(a;2, 3:3, &2, &3, "^exc) 
-X2/iexc(a;3,a;2,&3,&2,"^oxc)] , (3.13) 

with the constants Cint = Cexc = -Cint = -Cxc = 1- 

In Eqs. ( p.9|) , ( p.lO|) , ( p.l2|) , and ( p.l3|) the Wilson coefficients ai and 02 are defined by 



(3.12) 



X3 



ai = Ci 



9i 



^2 = C2 



9i 



(3.14) 
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with Nc = 3 the number of colors. The functions /I's, obtained from the Fourier transfor- 
mation of the lowest-order H, are given by 

h{xi,X2, &i, &2, rn) = Kq (y/x^x^bi) 

X [9{bi - 62)^0 iVx^bi) Iq {y/x2mb2) 
+e{b2 - b,)Ko (v/Si^&2) /o (V^i^&i)] , (3.15) 

hint{xi,X3, 61, 63, mint) = h{xi, Xs, bi, 63, mint) , (3.16) 

/iexc(a:2,a;3,&2,&3,"^cxc) = x^o^^ {Vx^x^rn^b2) 



X 



0{b2 - &3)-f^o (V^3"^cxc&2) ^0 i^xsnic^cbs 
9{bs - b2)H^^^ (y/x^m^bs) Jq {y/x^m^M) 



(3.17) 



with m = M^ and mint = "Tiexc = (1 ^ r^)M^. The hard scales t are chosen as 



t = max (y^a^Lm, y^aJim, 1/61, 1/62) , (3.18) 

tint = max(Va:imint, ^a^s^^int, I/&1, 1/&3) , (3.19) 

texc = max(^X2mexc, A/xirrw, I/&2, 1/&3) • (3.20) 

For the nonfactorizable amplitudes, the factorization formulas involve the kinematic 
variables of all the three mesons, and the Sudakov exponent is given hj S = Sb + Sjj(*) + Sj,-. 
The integration over 63 can be performed trivially, leading to 63 = 61 or 63 = 62- Their 
expressions are 



1 




M'CI = 327iV2NCfV^MI f [dx] rb^db^b2db2<PB{xiM)(pDi*){x2,b2)<PAxz) 



C (t ) 



X I [2:3(1 - r') - xi - Ct^X2{r - r2)]/i«(x„ 6,) 

-[^3(1 - r^) -x,+ X2]hf\x„ 6,)| , (3.21) 

Mtl = S2ttV2NCfV^MI I [dx] rbidbib2db2(t)B{xiM)(t)Di*){x2M)Mx^) 

Jo Jo 



xasitd) ^^'^ exp[-g(trf)|b3=fcj 
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x| [xi -X2 -X3(l - r'^)]h^J\xi,bi) 

-[(xi + X2)(l + dV) - l]hf\x., 6.)} , (3.22) 

A^W = 327rv/2iVCFV^M| / [dx] / 6it/6i62C?&20i3(a;i, 6i)0^(*)(x2, &2)0.(a;3) 



X"«(^/)^^exp[-5(t/)|b3=b2j 



A^ 



x}^[xs{l-r')-C^*\x^-X2y]hf\x„b,) 
-[(xi + X2)(l + Cfr^) - c}*V2]/if (x„ 6,) 

with the definition [dx] = dxidx2dx3. The constants are (b,d. 



(3.23) 



/ 



-c 



f>,d,/ 



1. 



The functions h^^\ j = I and 2, appearing in Eqs. (|3.21|) - (|3.23|) , are written as 

hi'^ = [e{h - h)K, {BMsh) Jo {BMBh) 

+eib2 - b,)Ko {BMBb2) /o {BMBbi)] 



Ko{B,MBb2) 



Aj) 



for 5| > 
\B[\MBb2) for 5| < 
[e{b, - b2)Ko {DMBbi) h {DMBb2) 



ITT 0-(l)l 



+9{b2 - b,)Ko {DMBb2) h (DMBb,)] 



X 



Ko{D,MBb2) 



for Dj > 



M) 



VT r. 



D]\MBb2) for D| < 

r(l) 



^(6i - b2)H'o'' (FMBbi) Jo (FMb&2) 
+0(62 - b,)Hl,'^ {FMBb2) Jo (FMb6i) 



KoiFjMBb, 



for F/ > 



r(i), 



\^fi7^^^(y'|i^2|j\^B6i) fori^2<o 



with the variables 



B = X1X2 , 

B"^ = xia;2 - a;2a;3(l - r^) 



Bi = XiX2(l + r^) - (xi - X2)(l - X3)(l - r^) 
D2 = a;,a;3(l_r2), 



(3.24) 



(3.25) 



(3.26) 
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DI = {xi + X2)r^ - (1 - xi - X2)x3(l - r^) , 

F' = X2Xs{l - r^) , 

F| = xi + X2 + (1 - xi - X2)x3(l - r^) . (3.27) 



The scales t^^^ are chosen as 



U = max{BMB, J\Bf\MB, J\BI\Mb, 1/&i, I/62) 



U = maxiDMs, J\DI\Mb, J\DI\Mb, I/61, 1/52 



tf = max(FMB, ^J\F?\Mb, ^\Fi\MB, l/h, l/h) . (3.28) 

As before, we have included the intrinsic b dependence of wave functions for the heavy 
B and D^*^ mesons, in which soft dynamics is more important. This dependence for the 
fast recoiling pion is also neglected, similar to the case of the K* meson, for consistency. 
The scheme of a single hard scale t for all diagrams contributing to an amplitude has been 
adopted to simplify the factorization formulas. We have also corrected misprints in Eq. (p. 27 ) 
appearing in 0. 

IV. NUMERICAL ANALYSIS 

In this section we determine the heavy-meson wave functions by fitting the predictions 
from the factorization formulas for the B — > K*'j and B -^ D^*^7r decays to the CLEO data 



B{B -^ K*-f) = (4.2 ± 0.8 ± 0.6) x 10^^ , 

B(B- -^ D\-) = (5.0 ± 0.5 ± 0.2) x 10"^ , 

B{B^ -^ D+n-) = (3.1 ± 0.4 ± 0.2) x 10"^ , 
B{B- -^ D*^Tx~) = (4.34 ± 0.33 ± 0.34 ± 0.18) x 10"^ , 
B{B^ -^ D*+n-) = (2.81 ± 0.11 ± 0.21 ± 0.05) x 10"^ . (4.1) 
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We adopt Gp = 1.16639 x 10"^ GeV-^, the CKM matrix elements \Vcb\ = 0.040 |||T|] and 



\V^^\ = 0.974, the masses Mb = 5.28 GeV, Md = 1.87 GeV, and M^* = 2.01 GeV |2T|, and 
the 5° {B-) meson hfetime r^o = 1.53 (r^- = 1.68) ps p2|. 



As the transverse extent b approaches zero, (j){x, b) reduces to the standard parton model 
0(x), i.e., 0(x) = (f){x,b = 0). The wave functions 4>i{x), i = B, D, D*, K*, and vr, satisfy 
the normalization 

f^ 



(j)i(x)dx _ 

^^ ' 2^6 



(4.2) 



with the corresponding decay constants /j. We choose fs = 180 MeV (in the convention 
/^ = 132 MeV), which is reasonable according to the lattice calculations presented in the 
literature [^]. For the wave function of the K* meson with transverse polarizations, we 
employ [110] 



\/6 

^K'{x) = ^fK'x{l -x)[l + 0.6(1 - 2x) + 0.06(5(1 - 2xf - 1)] , 

with /x* = 185 MeV. As to the pion wave function, we employ the asymptotic form, 

07r(a;) = — /^x(l - x) , 



(4.3) 



(4.4) 



which has been found to give predictions consistent with the data of the 717 form factor [24 



For the B meson wave functions, we use the model P5 



0b(x, 6) = Nb\J x{l — x) exp 



1 /xMb 

2 V ujb 



' ^|&^ 



(4.5) 



where Nb is the normalization constant and ujb the shape parameter. Another model with 
a fiat profile [||] 



(/)'s(x,6) 



N'Bbx\l xf ._KAJMlx + CBx{l-x)b 



(4.6) 



'Mix + Cbx{1 - x) 

gives predictions for the branching ratio of the B —>■ K*^ decay, which are too small com- 
pared with the data. Therefore, we shall not consider it in the numerical analysis below. For 
the D^*^ meson wave functions, we adopt the same model as of the B meson in Eq. ( [4.5[ ), 
but with different normalization constants A'^q(,) and shape parameters Wbx*), 
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0£,(*)(x,6) = Njj(,)^x{l -a;)exp 



1 (xMoi*)V ^l^b' 



(4.7) 



2 y t^£)(.) / 2 

We shall determine 5 parameters NoifD), N^'ifD*), ^b, ^d, and ujd, from the data 



of the 5 processes in Eq. ( [4.1| ), under the constraints /£>(.) > /b, which is reasonable from 
the viewpoint of heavy quark symmetry. The fitting procedures are as follows. Start with 
the K* meson wave function in Eq. ( [4.3|) , and determine the B meson wave function {i.e., 
the shape parameter ujb) from the data of the decay B -^ K*^. The obtained B meson 
wave function is then employed to evaluate the branching ratios of the charged and neutral 
B decays B -^ 0^*^71, from which both parameters, Njj{*) and u;£)(*), of the D^*^ meson 
wave function can be extracted. We have tested the heavy meson wave functions without 
the intrinsic b dependence. It turns out that the D'-*^ meson wave function 0£,(,)(x), along 
with the B meson wave function (J)b{x) determined from the B -^ K*'y data, always lead to 
predictions for the B -^ D^*^it decays which exceed the data by a factor 3, no matter how 
0£)(.)(x) is adjusted under the constraint f^x*) > /b- These predictions are lowered only by 
introducing the intrinsic b dependence into heavy meson wave functions, implying that the 
intrinsic b dependence is indeed essential. 

It is easy to understand that the parameters uj^(t) control the difference between the 
charged and neutral B meson decay rates. This fact is obvious from Eqs. ( p.5|) -( p^ : by vary- 
ing co'^(*), the external Pi^-emission contributions, i.e., the B ^ D form factors change, while 
the internal Pi^-emission contributions, i.e., the B -^ it form factors remain the same. The 
VT-exchange and nonfactorizable diagrams give smaller contributions than the VT-emission 
diagrams to the decays considered. Therefore, the charged B decay rates, as a combination 
of the external and internal Pi^-emission contributions, and the neutral S-meson decay rates 
with the external W emissions as the main contribution, vary differently with the change of 

UJb,(*). 



After a tedious numerical analysis, we determine the meson wave functions, 
0b(x, b) = l.U605Jx{l - x) exp -- ( ^ ^ \ - -(0.462 GeV)^6' 



2 V0.462 GeV/ 2 



0.3 



0.2 



0.1 



0.0 



1 1 






i\ 

1 \ 






C^o(x) 

^o.(x) 


9k.W 


r \ 




- 



0.0 



0.5 

X 



1.0 



FIG. 1. The K*, B and D^*' meson wave functions. 



0b(x, h) = 0.218677 Ja;(l - a;) exp 



1 f xMp Y 

2 V0.8 GeVy 



21,2 



-(0.8 GeV)'6' 



^D* {x, b) = 0.202997a/x(1 - x) exp 



Kofe)->-v)^^^' 



(4.8) 
= 190 



which correspond to the decay constants fs = 180 MeV, fo = 220 MeV, and fo* 
MeV, respectively. With these heavy meson wave functions, we obtain 

B{B -^ K*-f) = 4.19 X 10"^ 

B{B- -^ D\-) = 4.74 X 10"^ 

B{B° -^ D+TT-) = 3.50 X 10"^ 
B{B- -^ D*°n-) = 3.93 x 10"^ , 
B{B° -^ D*+TT-) = 2.94 X 10"^ , (4.9) 

which are consistent witli the data in Eq. ( ^.1| ). Note that the value ujb = 0.462 GeV is in 
good agreement with ub = 0.4 GeV obtained in the oscillating quark model |2^. It indicates 
that meson wave functions in PQCD factorization theorems have absorbed nonperturbative 
dynamics in B meson decays appropriately, and that remained dynamics, collected by the 
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hard quark-level amplitudes, is perturbative. The x dependence of the extracted heavy- 
meson wave functions and of the K* meson wave function is displayed in Fig. 1. As expected, 
the maximum of a heavy meson wave function moves to the large x region as the meson 
mass decreases. 

V. CONCLUSION 

In this paper we have extracted heavy meson wave functions, as shown in Eq. (|4.8| ), from 
the experimental data of the B — > K*'y and B -^ D^*^n decays. The extracted meson decay 
constants fo and fo*, compared to fs, and the behavior of the wave functions exhibited in 
Fig. 1 are all reasonable. We have also confirmed that the intrinsic b dependence of heavy 
meson wave functions, which makes possible a simultaneous fit to the data of several decay 
modes, is essential. Although the uncertainties of the data are still large, such that wide 
ranges of the relevant parameters are allowed, we shall employ the obtained wave functions 
to make predictions for other nonleptonic B meson decay modes, especially charmless de- 
cays. As the precision of data improves, the heavy meson wave functions determined in 
our formalism will be more reliable. On the other hand, we can also proceed with a global 
determination of the heavy meson wave functions by incorporating more decay modes into 
the fitting procedures. These subjects will be discussed elsewhere. 
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